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Theorems of Euclidean Geometry through Calculus 
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We re-derive Thales, Pythagoras, Apollonius, Stewart, Heron, al Kashi, de Gua, Terquem, Ptolemy, 
Brahmagupta and Euler’s theorems as well as the inscribed angle theorem, the law of sines, the 
circumradius, inradius and some angle bisector formulae, by assuming the existence of an unknown 
relation between the geometric quantities at stake, observing how the relation behaves under small 
deviations of those quantities, and naturally establishing differential equations that we integrate out. 
Applying the general solution to some specific situation gives a particular solution corresponding to 
the expected theorem. We also establish an equivalence between a polynomial equation and a set 
of partial differential equations. We finally comment on a differential equation which arises after a 
small scale transformation and should concern all relations between metric quantities. 


I. THALES OF MILETUS 


Imagine that Newton was born before Thales. When 
considering a triangle with two sides of lengths x and 
y, he could have fantasized about moving the third 
side parallel to itself and thought: ”Well, I am not 
an ancient Greek geometer but I am rather good in 
calculus and I feel there might be some connection 
between the way x and y vary in such circumstances.” 


He would have materialized 
his suspicion in a function 


y = y(z) (1) 


connecting x and y what- 
ever the position of the third 
side, as long as it is moved 
parallel to itself. In partic- 
ular, after a slight displace- 
ment resulting in small devi- 
ations 6x and dy, he would 
have had at first order 


dy=y'(x)dx (2) 


But the lengths da and dy 
of the small added segments 
must themselves obey equa- 


tion (1), that is 


dy = y(dx) (3) 


6x oy 


6x oy 


To see it, translate those segments to the (x,y) vertex. 
Developing the right-hand-side member of equation 
at first order and noticing that y(0) = 0, we have 

dy = y' (0)da (4) 
which, compared to eq. (2), implies that y’(x) is con- 
stant. Integrating y’ = k, k being a positive constant 


since y(a) is an increasing and smooth function, gives 
the Thales theorem [1] Book VI, Prop. I]] 


y(a) = ka (5) 
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II. PYTHAGORAS OF SAMOS 


If he was born before Thales, Newton was born before 
Pythagoras too, so that we do not have to make any 
further unlikely hypothesis. Imagine that driven by his 
success in suspecting the existence of a Greek theorem, 
he moved to consider a right triangle of legs of lengths x 
and y and of hypotenuse of length z. 


He might have been tempted 
to speculate about the link, 
if any, between x, y and z 


in every right triangle. And y 7 
again, as calculus master, he 
could have postulated that 
xz 
z= 2(x,y) (6) 


a relation that must be true for any xz, y and z ina 
right triangle. In particular, after a slight increase in the 
length of x, leading to a small deviation 6x, while dy = 0, 
he would have found, at first order, that 


bz = 0,26x (7) 


6z 
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Here, dx is the length of the hypotenuse of a right trian- 
gle with one leg of length 6z. At first order, this small 
triangle is similar to the initial one. Using the Thales 
theorem that he has just found out, 


bz = 6a (8) 
z 


Substituting this result to dz in eq. (7) would have led 
him to the partial differential equation 


20,2 = 2 (9) 
whose general solution is 


2 (x,y) = 2° + k(y) (10) 


where k(y) is an arbitrary function of y. But the function 
z has to be symmetric in z and y (i.e. he could have made 
the same reasoning with a non-zero dy while 6x = 0) and 
z(x,0) =a. Hence 


Aaya +y (11) 


that is, the Pythagorean theorem, which, like Thales’, 
was probably discovered long before — and published a 
few centuries later by Euclid [I] Book I, Prop. XLVI]. 
As will be shown in the last section, the same result can 
be obtained by considering a small rotation of the hy- 
potenuse around one of its extremities followed by an 
infintesimal scale transformation. This proof is known 
and is published in a slightly different form in [2] 3). 


Ill. APOLLONIUS OF PERGA 


Thales and Pythagoras the- 
orems are not the only 
ones that are named _ be- 
fore famous Greek geome- 
ters. Newton could have 
gone a step further — east- 
wards, a few centuries later 
— and assumed that in any 
triangle of sidelengths x, y 
and z, the length d of the 
median relative to the z- 
length side is a smooth func- 
tion of x, y and z, ie. 


d= d(x, y, z) (12) 


After an infinitesimal rota- 
tion of the y-length side 
around the (y,z) vertex re- 
sulting in a small deviation 
da, with dy = 0 and dz = 0, 
at first order: 


dd=0,d5x (13) 


In order to get an expression for dx and 6d and then 
a differential equation leading to the would-be theorem, 
consider the infinitesimal arc travelled by the moved ver- 
tex, of length d¢. 


At first order, it can be seen 
as the hypotenuse of a small 
right triangle with one leg 
of length 6x, which is sim- 
ilar to a larger right triangle 
whose corresponding leg is 
the hy-length height relative 
to the y-length side, and the 
hypotenuse is the z-length 
side, so that 


The délength arc is also 
the first order hypotenuse of 
another small right triangle 
with one leg of length dd, 
which is similar to the tri- 
angle with a d-length hy- 
potenuse and whose corre- 
sponding leg is a segment 
starting from the foot of the 
median and parallel to — and 
thus half of the length of — 
the hy-length height, so that 


hy 
bd = S4oe (15) 


Inserting those deviations in eq. leads to the partial 
differential equation 


dd,d = 5 (16) 


whose general solution is 


P(x,y,2) => +hY,2) (17) 


where k(y, z) is a function of y and z. But d(x, y, z) has 
to be symmetric in x and y (i.e. we can make the same 
reasoning with a non-zero dy while dx = 0). Hence 


2 2 
@(x,y,2) == 24 


+ ¢(z) (18) 


with c(z) a function of z. Furthermore, if « = 0 (or 
y =0), y=z (or e = z) and d= 2/2. This yields c(z) = 
—z?/4, which can alternatively be found by invoking the 
Pythagorean theorem for y = x. The particular solution 
reads 


x? + y* — 2(z/2)? 
2 


that is, Apollonius’s theorem, to be found in a slightly 
more elaborate form in [4]. 


dé’ (a, ¥, 2) = 


(19) 


IV. MATTHEW STEWART 


Suppose Newton was born before Stewart, an 18th- 
century Scottish mathematician (and reverend). Well, 
he was. Perhaps he was not interested, or did not 
have the time, otherwise he could have used this 
tool to generalize Apollonius’s theorem to any cevian. 


In a triangle of sidelengths x, y and z, assume that the 
length d of a cevian dividing the side of length z in two 
segments of lengths m and n, is a smooth function of z, 
y, mand n, that is 


d= d(x,y,m,n) (20) 
After an infinitesimal rotation of the y-length side around 
the (y,z) vertex resulting in a small deviation 62, with 


dy = 0, 6m = 0 and 6n = 0, at first order: 


5d = 0,d 6x (21) 


Using the same similarities as for the Apollonius’s theo- 
rem, with the unique difference that the foot of the cevian 
is not necessarily the middle of the (m+ n)-length side 
but falls at a distance n from its right vertex, we find 


ery ae eee 


x (m+ aa 


hy being the length of the height relative to y, we have 


the partial differential equation 


n 


ddzd = 


x (23) 


m+n 


whose general solution is 


n 


d*(x,y,m,n) = x? +k(y,m,n) (24) 


men 


where k(y,m,n) is a function of y, m and n. But 
d(x,y,m,n) must be symmetric in (x,m) and (y, 7) (ie. 
we can make the same reasoning with a non-zero dy while 
dx = 0). Hence 


nx? + my? 


d?(x,y,m,n) = + c(m,n) (25) 


m+n 


with c(m,n) a symmetric function of m and n. Further- 
more, if « = 0 (or y = 0), y= m+n (or x = m+n) 
and d= m (or d=n). This yields k(m,n) = —mn. The 
particular solution reads 


n(a? — m?) + m(y? — n?) 


m+n 


d(x,y,m,n) = (26) 


that is, Stewart’s theorem [5]. 


V. HERON OF ALEXANDRIA 


Intoxicated by his findings, Newton could have switched 
to a more elaborate, though older, challenge — as proba- 
bly did an Ancient Greek Roman Egyptian mathemati- 
cian... What if, for any triangle, the area A could be 
a smooth function of the sides lengths x, y and z? He 
would have assumed 


A= A(z, y,z) (27) 


After an infinitesimal rotation of the y-length side around 
the (y, z) vertex resulting in a small deviation 6x, while 
dy = 0 and 6z = 0, at first order: 


5A =0,A6x (28) 


First note that the h,-length 
height relative to the z- 
length side divides the ini- 
tial triangle in two right 
triangles of horizontal legs 
of lengths ¢ and z —t re 
spectively. One can ex- 
press hz, as a result of 
the Pythagorean theorem in 
both right triangles. 


Equating those expressions yields x? — t? = y? — (z —t)? 
and hence 
wo? — y? + 22 =a ee? 
t= z-t= 
22 22 


(29) 


bx_ 6é 


toh pz-t 
Zz 


z 
Again, 0@ is the length of the infinitesimal arc travelled 
by the moved vertex. Like in the two last sections, it can 
be considered as the first-order hypotenuse of a small 
triangle whose similarity with a larger one allows to find 
6a. But it is also, at first order, the hypotenuse of another 
small triangle with one leg of length dh,, which is similar 
to the large right triangle whose corresponding leg is the 
(z — t)-length segment, and the hypotenuse the y-length 
side. Since hy = 2A/y and 6h, = 26A/z, we have 


y? — 72 4 72 
dy 


Plugging in results into equation (28), gives the par- 
tial differential equation 


ee aS 5€ (30) 
zy 


Ad, A = = [o(y? + 22) — 25) (31) 


which can be integrated out to give the general solution 


1 
A*(x, 9,2) = Fe [2e*(y? + 2°) — x +k(y,z)| (32) 
where k(y,z) is an homogeneous function of y and z. 
Since A(x, y, z) must be symmetric in x, y and z (i.e. we 
can make the same reasoning with a non-zero dy or 62), 
k(y, z) = 2y2z? — y* — 2+. Hence 


1 
A(t, 942) = GVXOP FPP EPP) — yee) 
(33) 


which can be factorized into the Heron theorem [6] 
ytezu+y-—Zz 
2 2 


e+ytz-H+y+z2 
A(z.y.2) = yf 5 5 
(34) 


whose discovery could actually be Archimedes’ [7]. 


VI. JAMSHID AL-KASHI 


Newton could have chosen to deal with angles — besides 
calculus, he knew a bit about trigonometry. Let us send 
him to Persia, a few centuries before his birth, and won- 
der wether in any triangle of sidelengths x, y and z, the 


_—__ 


angle 7 = (x,y) could be a smooth function of x, y and 
z, that is 


Y= (2,92) (35) 


After an infinitesimal rotation of the y-length side around 
the (y,x) vertex resulting in a small deviation 6z, with 
dx = 0 and dy = 0, at first order: 


Oy = O27 62 (36) 
6z_ 5é 
z y z 
oy 
x Y x 


While dy is easy to connect to 6@, the length of the arc 
travelled by the moved vertex (in the illustrative figures, 
x and z have been swapped for aesthetic reasons), dz can 
be determined thanks to the same similiarity as in the 
three previous sections. We have thus 


Mt 
y 


Inserting those deviations in eq. yields the partial 
differential equation 


h 
oy 6z= 6 with hy=azsiny (37) 
z 


siny 0.7 = = (38) 
ty 


whose general solution is 


27° +k(a,y) 
ee aa 39 
cosly(, 2) = (39) 
where k(x, y) is a symmetric, homogeneous function of x 
and y. According to Pythagoras, when y = 7/2, 27 = 
x? +y?, ie. k(x, y) = —x? — y?. Hence 
2 4er+y? 
= 40 
cosly(e,y,2)] = = (10) 
that is, al-Kashi’s theorem [8] — also known as the law of 
cosines or generalized Pythagorean theorem, and already 
familiar to Euclid [I] Book II, Prop. XH & XI]. 


VII. OLRY TERQUEM 


Completely exhilarated, 
Newton could have taken on 
a bigger piece and assumed 
that in any triangle of 
sidelengths x, y and z, the 


: \ length d of the y = (2,y) 
angle bisector is a smooth 

a function of x, y and z, i.e. 
; SN, dada.) (4D 


After an infinitesimal rotation of the y-length side around 
the (x,y) vertex resulting in a small deviation 6z, with 
da = 0 and dy = 0, at first order: 


bd = 0.d6z (42) 


Again, thanks to the same similiarity as in the four pre- 
vious sections, 6z can easily be linked to 62@, the length 
of the arc travelled by the moved vertex. 


It is a little more complicated for dd. First note that in 
the illustrative figure, dd < 0, so that we will consider 
the positive length —dd. Then observe that when the y- 
length side infinitesimally rotates around the (a, y) ver- 


tex, the foot of the y = (a, y) angle bisector moves along 
a perpendicular to the y-length side, just like the (z, y) 
vertex. But the angle between this perpendicular and 
the angle bisector is the complementary of y/2. Thus in 
the small right triangle of legs of lengths —dd and ddéy/2, 
the opposite angle to the —dd-length leg is, at first order, 
equal to 7/2, implying that tan(y/2) = —dd/(ddy/2). 
Hence 


d 
oz = 00 dd = — tan 457 with dé = ydy (43) 


z 


hy being the length of the height relative to y. Using 


tan 2 = 927 _ (44) 
2 1+ cosy 
with 
hy _ 42 2442 
siny=—% and cosy= — os (45) 
x 2ry 
we have the partial differential equation 
Od — 
= (46) 
d~ =2+ +p 
whose general solution is 
d(x, y,z) = k(x, y)/(@+y)? — 2? (47) 


where k(x,y) is a symmetric function of 2 and y. To 
determine it, note that in the particular case of a right 
triangle with hypotenuse of length z, the angle bisector is 
the diagonal of the inscribed square of sidelength xy/(a+ 
y) — as can be deduced from similarities between the right 
triangles generated by the square in the initial triangle. 
We find k(x, y) = ,/xy/(x + y). Hence 


d(x, y,2) = = (1 = oa) (48) 


that is, the length of the angle bisector, as Terquem com- 
puted in the 19th century [9]. 


VIII. JEAN-PAUL DE GUA DE MALVES 


Armed with this powerful theorem-finding tool, Newton 
could have moved on to even bolder challenges, like leav- 
ing the plane for the real space, and imagining, say, a 
generalization of the Pythagorean theorem in three di- 
mensions! Let him consider a trirectangular tetrahedron, 
that is a tetrahedron with a right angle corner, like the 
corner of a cube: what if, for any of them, the area of 
the face opposite to the right angle was a function of the 
areas of the other faces? 


5 


A convenient way to 
parametrize the problem is 
to give arbitrary lengths to 
the three edges from the 
right angle vertex, say 2, 
y and z. The areas of the 
three right triangle faces 
are ry/2, 22/2 and yz/2. 
For the area of the last 
face, opposite to the right 
angle, say A, we can have 
an expression by choosing a 
base, say the edge of length 
/y? + 22 (thanks Pythago- 
ras) and the relative height 
of length h. We have 


A=5VP Teh (49) 


Let us go back to Newton and his obsession. He could 
have stated that A is a smooth function of x and y: 


A= A(z, y, z) (50) 


Choosing to slightly increase x, while leaving y and z 
invariants, that is, an infinitesimal deviation 62, with 
dy = 0 and 6z = 0, we find 


=o Ase (51) 
Eq. (49) implies that 


bA= VP +e oh (52) 


But what do we know of dh? First note that the foot 
of the h-length height is not affected by the deviation 6x 
since this h-length height and the z-length edge are in a 
plane orthogonal to the base of the A-area face. In this 
plane, we can check that at first order, the right triangle 
with h-length hypotenuse and z-length leg is similar to 
the one with dx-length hypotenuse and dh-length leg, so 
that 


bh = bx (53) 


Combining this equation with result (52), itself plugged 
in into eq. with dy = 0, we have 


VP +2 be = 0,A 6x (54) 


Simplifying by dx and using eq. (49) to get rid of h, we 
find a partial differential equation 


AQ, A = ACs Phe (55) 


It can be integrated out to give the general solution 


A’(az,y,2) = Gly? + 27)2? + k(y, 2)] (56) 


1 
4 
with k(y, z) an homogeneous and symmetric function of 
y and z. Since A(x,y,z) must itself be symmetric in a, 
y and z (ie. we can make the same reasoning with a 
non-zero dy or 6z), k(y,z) = y2z7. Hence 


sons) (2) +B) +(B) on 


known as de Gua’s theorem [10], first formulated by 
Descartes [11], which states that in any trirectangular 
tetrahedron, the square of the area of the face opposite 
to the right corner is equal to the sum of the squares of 
the areas of the other faces — a three-dimensional gener- 
alization of the Pythagorean theorem. 


IX. THE INSCRIBED ANGLE 


Let us move on to the circle, and confront Newton to a 
simple problem. 


Assume that an arbitrary 
inscribed angle a is a 
smooth function of the cen- 
tral angle @ that intercepts 
the same arc on the circle, 
ie. 


) a=a(0) (58) 


After an infinitesimal devi- 
ation 60, at first order, 


da=a'd™ (59) 


At first order again, the 
small right triangle with 
one leg of length yéa and 
hypotenuse of length 62 is 
similar to the triangle of 
corresponding sidelengths 
y/2 and R respectively. 


Thus we have 


you = —— de (60) 


Since 60 = 6€/R, this leads to the differential equation 
a’ = = (61) 


whose general solution is 


a(0)=S +k (62) 


where & is a constant that vanishes since a(0) = 0. Hence 
a(?) = 5 (63) 


that is, the inscribed angle theorem. 


X. THE CIRCUMRADIUS 


We could keep Newton in the circle and think to another 
question: of course, for any triangle, the circumradius 
length R should be determined by the sidelengths x, y 
and z. It was known to Euclid {1} Book IV, Prop. V] and 
it is actually simple to prove it with the help of Heron’s 
theorem. But let him play the game of finding it from 
scratch, that is, by postulating that 


R= R(a,y, z) (64) 


£ 


Li 


If the (x, y) vertex is slightly moved along the circumcir- 
cle, it generates infinitesimal deviations 6x and dy, while 
6z =O and 6R=0. At first order, we have: 


0,R dx + OyRdy = 0 (65) 


We need to go through 
some geometric consider- 
ations before proceeding: 
any angle in the triangle 
has the same magnitude as 
the ones between its op- 
posite side and the tan- 
gent lines to the circumcir- 
cle from the two other ver- 
tices, since they intercept 
the same arc. 


After the infinitesimal displacement of our vertex along 
the circumcircle, we consider the right triangles with dé 
length hypotenuse and dz-length leg, and with 6¢-length 
hypotenuse and —dy-length leg, respectively, and derive 


da = cosa dé — dy = cos Bd (66) 


Applying the al-Kashi theorem yields 


2 2 2 2 2 2 
YE ty +e 
= as Je iy=—z (67) 


Plugging in results into eq. (65), simplifying by 
6€/2z and isolating O,R gives the partial differential 
equation 


6x 


2 A242 
vy —-xw+z 
A similar expression arises for 0, R when the (, z) vertex 
is slightly moved along the circumcircle: 


2 2 2 
=, + 
LZ x Yy R 


0,R = -— >—— 5 
@ za2—224y? 


(69) 


Finally, we consider the infinitesimal scale transformation 
Rt R+6R. Since each of the sidelengths increases in 
line with itself and 6R/R, that is dx = «/ROR, etc., we 
have the partial differential equation 


KR=20,8+y0,8-+ 20,8 (70) 


which guarantees that each quantity at stake is taken 
into account with its appropriate dimension. Plugging in 


eq. and into the last one yields the partial 


differential equation 


OrR g* — y* — 24 + Qy? 2? 
R a[—a4 — yt — 24 + 2(y22? + 2222 + a2y?)] 


(71) 


which can be integrated out to give the general solution 


xt k(y, 2) 
R(x, ¥, 2) = 4 4 4 252 2 52 27,2 
J—a! — yt — 24 + 2(y2z? + 222? + xy?) 
(72) 
with k(y,z) a symmetric function of y and z. Since 


R(x, y, z) must itself be symmetric in x, y and z, k(y, z) = 
cyz where c is a constant. We just have to compute 
the circumradius length of, say, an equilateral triangle 
of sidelength 1, that is R(1,1,1) = V3/3, to find out 
that c = 1. Furthermore, like for Heron’s theorem, the 
argument of the square root can be factorized to give 


R(x, y, 2) (73) 
LYz 
V(at+yt+z)\(—a@t+yt2y(x@-yt2z\(@t+y-2) 


which is indeed the expression of the circumradius as a 
function of the sidelengths. 


XI. THE LAW OF SINES 


Let us go back to the Thales theorem that Newton could 
have discovered by the grace of calculus if only he was 
born two millenia before: y(x) = kx is already a nice 
result, that we have used all along this paper, but what 


about k? If he knew the angle y = (x,y) between the 


a 


two sides, and one of the angles, say 6 = (z,x), between 
the z-length side and the third side, of length z, he could 
have speculated that y is a smooth function of x, 6 and 
y, that is 


y = y(x, B,) (74) 
It is possible to find the 
function y with a sim- A 


ple rotation of one side 
around one of the adja- 
cent vertices, like in sec- 
tions [ITI] to but the © 
differential equations that 
arise are a bit complicated 
to solve: one is linear 
and the other one of the 
Bernouilli type. Even if 
it is not straightforward, it 
happens to be much sim- 
pler to work within the cir- 
cumcircle. After an in- 
finitesimal and clockwise 
displacement of the y ver- 
tex along the circumcircle, 
resulting in a small devia- 
tions dy, da and 68 while 
oy = 0, at first order, 


dy = Ozy Ou + Ogy 68 (75) 


The deviations 6a and dy are the same as in section [X] eq. 
(66). —08 can be computed by considering —27d6 as the 
length of the leg opposite to a in the small right triangle 
of other leg of length 6a and hypotenuse of length 6. 
Using a = 7 — (8 +7) to get rid of a, we have at first 
order 


_sin(B +7) sp 


(76) 
Combining eq. and (76) leads to the differential 
equation 


—dy=cos8dl dx=—cos(8+y7)de —68 


xcosB = xcos(h+ 7)0zy + sin(B + y)Ogy (77) 


If we now consider an infinitesimal and clockwise dis- 
placement of the 6 vertex along the circumcircle, result- 
ing in deviations dx and dy while dy = 0 and dy = 0, at 
first order 

O,y 6x + O,y dy = 0 (78) 
Using the analogous first-order similarities and the fact 
that a = 7 — (8+-7), we find 


— dx = —cos(8 + 7) 6 oy = mee) 7) b€ (79) 
Combining eq. and (79) leads to the differential 


equation 
xcos(8+7)dzy + sin(6 + 7)d,y = 0 (80) 


Finally, we consider the infinitesimal scale transformation 
R+> R+06R, where R is the circumradius. Since each 
of the sidelengths increases in line with itself and 6R/R, 
that is dx = «/ROR, but 68 = 0 and dy = 0, we have 
the partial differential equation 


y = xOgYy (81) 


Using eq. (81) allows to simplify eq. and to 
yield 


y _ «cos B 

Oey tan(8+y)  sin(8 +7) 482) 
y = 

caer a a on 


They are both linear. Moreover, the second one is sepa- 
rable and admits the general solution 
K (a, B) 

7 7 — . 84 

y(x, B,) sin(8 +7) (84) 

where AK (a, 8) is a function of x and 8. But this solution 

must satisfy eq. (81), that is K(«,8) = k(8)a where 

k(8) is a function of 8. The solution must also satisfy 

eq. (82), thus we find k(8) = sin8+C with C a real 

constant. Finally, since for 6 = 0, y = 0, we have 


xsin Gb 
sin(6 + 7) 


which, using sin(G+) = sina, is the law of sines, proved 
by Nasir al-Din al-Tusi, Persian mathematician — and 
architect — regarded as a founder of trigonometry [14]. 


y(z, 8, y) = (85) 


XII. CLAUDIUS PTOLEMY 


We now take Newton back 
to Antiquity, in Roman 
Egypt, and make him look 
at an ancient problem: 
given a cyclic quadrilat- 
eral, determine the length 
of its diagonals as func- 
tions of the sidelengths z, 
y, u and v, that is, for the 
z-length segments joining 
the (v,xz) and (y,u) ver- 
tices: 
a= lL, Y; U, v) 


After a small displacement of the (x, y) vertex along the 
circumcircle, resulting in deviations dx and dy while du = 
dv = 0 as well as dz = 0, we have at first order 


0,2 6% + Oyz dy = 0 (87) 
The deviations 6x, dy are those of section [X| eq. (66): 


dx = cosa dl — dy = cos Bd (88) 


where a and @ are the angles opposite to the x and y- 
length sides of the (x,y, z) triangle. 


Plugging in those results into eq. (87) and considering 
the equivalent after a small displacement of the (u,v) 
vertex, we find the two partial differential equations 


cos a 0zz — cos 8 Oyz =0 (89) 
cos Pp Oyz — cos Oyz = 0 (90) 
> where p and @ are the 


angles opposite to the u 
and v-length sides of the 
(u,v, Z) triangle. 

If we now slightly move 
the (y, wu) vertex along the 
circumcircle, with small 
deviations dy, du and 62, 
while dx = dv = 0, we have 
at first order 


ds 


6z = Oyz dyt Oz du (91) 


The deviations dy, du and 
dz are given by 


dy = cos Bd 
—du = cos pol 
—dz = cos yoel (92) 


where dy and —du are 
calculated analogously to 
dx and —dy above, while 
—dz is to be seen as the 
length of one of the legs 
in a small right triangle 
of hypotenuse of length 6@ 
and adjacent angle of same 
magnitude ¥ as the one be- 
tween the x and y-length 
sides. 


Inserting those results into eq. yields 
cos y = cos p0,z — cos 6 Oyz (93) 
Moreover, from the infinitesimal scale transformation 
Rw R+O0R, a fourth differential equation arises: 
Z=LOpz + YOyZ + UOyZ + V Oyz (94) 


Combining equations (89), (90), and (04), isolating 
O0,z and using «cos 8+y cosa = z and ucosp+vucos¢ = 
z to simplify the expression, we have 

(cos pcos @ — cos y) cos 3 


On2 = 95 
. cos pcos d + cos a cos 3 (5) 


Looking at the intercepted arcs, we observe that (u,v) 
and y are supplementary. We then exploit y = 7 — (a+ 
8) =p+¢ and the angle sum identities to find 


sin psin ¢ cos 3 


= 96 
sin psin@ + sina sin B 8) 


Thanks to the law of sines applied to both (2, y,z) and 
(u,v, z) triangles, and the al-Kashi theorem applied to 
cos 8, this equation finally becomes 

uv(x? — y? + 2”) 


(uv + vy)x 


220,2 = (97) 


and turns out to be linear in z?. Its general solution reads 
uv(a? + y?) + c(y, u, v)x 


(98) 


where the function c(y, u,v) must be symmetric in u and 

v, while z?(z, y, u,v) must notably be symmetric in x and 

y. Thus c(y,u,v) = y(u? + v). Hence, after factoriza- 

tion, 

(ua + vy)(va + uy) 
uv + ry 


27(2,Yy,U,v) = (99) 


which, if multiplied by the square of the other diagonal 
length, yields the theorem of Ptolemy [15]. Eq. (99) 
might have been first stated by Brahmagupta §28]. 


XII. BRAHMAGUPTA 


A few centuries later in India, was solved the problem 
of expressing the area A of any cyclic quadrilateral as a 
function of its sidelengths x, y, u and v. Well it could 
have been solved through calculus too. Let us start from 


A= A(x, y,u, 0) (100) 


Sw 


Like in the last section, we draw the z-length diagonal, 


which divides the quadrilateral in two triangles of sides 
of length x, y, z, and u, v, z respectively. Selecting the 
y-length side as the base of the (2, y, z) triangle, we com- 
pute its area by multiplying y and x sin y divided by two, 
where y¥ is the (a,y) angle magnitude. But the inter- 
cepted arcs of the circumcircle tell us that_in the other 


triangle, (u,v) is the supplementary of (x,y). Thus we 
have for their respective areas Azy, = xcysiny/2 and 
Aw: = uv siny/2. This yields 


2A 


——— 101 
ry + uv an 


siny = 
where A = Azyz + Auvz is the total area of the quadri- 
lateral. 

Let us go back to the project. Rather than moving 
the (a,y) vertex along the circumcirle, which leads to 
complicated equations, we invite Newton to follow the 
reasoning developed in section [V] for the Heron theorem 
and make an infinitesimal rotation of the y-length side 
around the (y,z) vertex resulting in a small deviation 
6x, while dy = 0 and dz = 0. We have at first order 


6A=0,A6x 


(102) 


This slight departure from 
the cyclical nature of the 
quadrilateral will be of no 
consequence. The reason- 
ing is valid in the triangle 
and as soon as the needed 
result is secured, we will 
get back into the circum- 
circle to proceed. Here h 
will stand for the length of 
the height relative to the z- 
length side in the (2, y, z) 
triangle, such that Azy, = 
zh/2 and thus 6A = zdh/2 
since dAyy2 = 0. We 
know a little more about 
trigonometry than in sec- 
tion [V] so that we can ex- 
press, at first order, the de- 
viations 6x and dh in terms 
of y and a, the (y, z) angle 
magnitude: 


dz = siny dé dh = cosa dé (103) 


where 6€ is the distance travelled by the moved vertex. 
Plugging in these into eq. (102), using eq. (101) for the 
sine and al-Kashi theorem for the cosine, we find 
Ab, A = TUT Pa 2 (104) 
8y 
that is, applying eq. to express z? in terms of the 
sidelengths, 


Ad, A = sea + (y? +u? + v7 )ax + 2yuv] (105) 


which can be integrated out to give the general solution 


A? (x,y, u,v) (106) 


= aa + 2(y? + u? + v?)x? + 8ryuv + ky, u, v)] 
where k(y, u,v) is a homogeneous function of y, u and v. 
A(x,y,U,v) must be symmetric in x, y, u and v. Conse- 
quently k(y, u,v) = —y*—u4—v4*#42(y2u? + y20? +-u?v?). 
After factorization, we come across the Brahmagupta 
theorem [16] §21] 


A(x, y, u,v) (107) 


_ (tee ty e—ytutv xty—utv r+ytu—v 
7 2 2 2 2 


XIV. LEONHARD EULER 


This one has a simple calculus part, but is a bit more 
trickier for the geometric principles. Imagine that New- 
ton, or even Euler, would have been interested in de- 
termining the distance d between the circumcenter and 
incenter of any triangle. Let him assume that d depends 
on the inradius r and the circumradius R, i.e. 


(108) 


After a slight displacement of one of the vertices, say 
O, along the circumcircle, resulting in a second triangle 
with small deviations dd and dr, while 6R = 0, we at first 
order 


bd = 0,d dr (109) 


Newton could not ignore that the circumcenter is the 
intersection of the perpendicular side bisectors and that 
the incenter is the intersection of the internal angle bisec- 
tors. Some geometrical considerations need to be taken 
into account in order to proceed. 


10 


First of all, the perpendicular side bisectors and the 
internal angle bisectors intersect on the circumcircle. To 
be convinced, just notice that the angles of the triangle 
are inscribed angles in the circumcircle. According to the 
inscribed angle theorem, their magnitude is determined 
by the intercepted arc. Hence, any angle bisector bisects 
the intercepted arc. So does the corresponding perpen- 
dicular side bisector. Let us call Q the intersection of 
the angle bisector from O and the perpendicular bisec- 
tor of its opposite side. Obviously, this angle bisector 
still passes through Q after the O vertex has been moved 
along the circle. 

Secondly, when a vertex is moved along the circumcir- 
cle, the incenter moves along another circle whose center 
is precisely the intersection of the angle bisector from 
the moved vertex, and the perpendicular bisector of the 
opposite side, that is, Q. 

It is simple to show that the trajectory of the incenter 
is circular: since the angle of the moved vertex remains 
constant throughout the movement, so does the angle of 
the moved vertex of a new triangle formed by the two 
remaining fixed vertices and the incenter; the inscribed 
angle theorem again tells us that the incenter moves along 
a circle passing through the two fixed vertices. 

To see that the center of this circle is on the circumcir- 
cle — and corresponds to Q — just move the incenter along 
its circular trajectory until it lays on the perpendicular 
bisector; in this symmetric configuration, a simple angle 
chase teaches us that the distance between the incenter 
and Q is the same as between each one of the fixed ver- 
tices and Q, so that this distance must be the radius of 
the circle, and @ its center. 

As a consequence, since it moves along a circle whose 
center is on the angle bisector from the moved vertex, at 
first order the incenter necessarily moves perpendicularly 
to this angle bisector. 

We can now go back to equation and try to find 
an expression for dd and dr by observing the way the 
initial triangle slightly moves to the second one. 


Note that dd (positive in the illustrative figure) is one 
of the legs of a right triangle whose hypotenuse length 
is the distance, say dp, between the incenters of the two 
triangles. Since at first order, the dp-length segment is 
perpendicular to the angle bisector from the moved ver- 
tex, and that the other leg is perpendicular to d-length 
segment, this right triangle is similar to another one of d- 
length hypotenuse and useful leg of length, say, g. Hence 


a 
dd = 7 op 


(110) 


For dr (negative in the illustrative figure), we need to 
draw the inradius perpendicularly to the side opposite to 
the moved vertex, and observe the way it evolves. Here, 
it is —dr that can be seen as the length of one of the 
legs of a right triangle with dp-length hypotenuse, which 
is similar to another one of hypotenuse of length R and 
corresponding leg of length g. We have 


br = — 5 6p (111) 


Dividing equations (110) and (111) member to member 


and using equation (109) yields the differential equation 


do,d = —R (112) 
which we integrate to obtain the general solution 
d?(r, R) = —2Rr + k(R) (113) 


with k(R) > 2Rr areal function of R. Since in an equilat- 
eral triangle, d?(R/2, R) = 0, we find k(R) = R?. Hence 


d?(r, R) = R(R—-2r) (114) 


that is, Euler’s theorem [12], first proved by Chapple [13]. 


11 
XV. THE ANGLE BISECTOR 


We already have a section dedicated to the angle 
bisector. But it was about its full length: from the vertex 
to the foot. What if Newton got now interested to the 
part between the vertex and the incenter? It should also 
be a smooth function of the sidelength x, y and z, i.e. 

c=c(z,y,z) (115) 
Since the position of the 
incenter matters, and that 
we know how it moves in 
the circumcircle, we will let 
Newton continue his con- 
siderations in this very cir- 
cumcircle. After a small 
clockwise displacement of 
the (z,y) vertex along the 
circumcircle, we have in- 
finitesimal deviations dy 
and 6z, while dz = 0 (in 
the illustrative figures, x 
and z have been swapped 
for aesthetic reasons). At 
first order, 


dc = Oycdy+0,cdz (116) 


In the small right triangle with legs of length —dc and 


coy/2 — where y = (x,y) — the angle opposite to the 
céy/2-length side has the same magnitude as the one 


—. 


between the z-length side and the 8 = (a, z) angle bi- 
sector, that is 6/2, since they both are inscribed angles 
intercepting the same arc of the circle along which the 
incenter travels. Thus we have 


B _ cby/2 
a 
2 —dc 


(117) 


Furthermore, based on the same reasoning as in the last 
section, we can connect —dy and 6z to d€, as well as 
yoy = sin G62, so that using eq. for the tangent, we 
find 


—c(1+ cos £) 


— dy=cosBdl dz=cosydl dbc= 5 
y 


be 


(118) 


Plugging in results into equation gives a par- 
tial differential equation whose analogue can be obtained 
by considering a clockwise rotation of the (x, z) vertex 
instead of the (z,y) one. Both equations read 


Cc 
2y 


51 + cosa) = cos a 0;c — cosy 0z¢€ 


(1+ cos 8) = cos 8 Oyc — cosy Oz¢ (119) 


(120) 


Combining those with the differential equation arising 
from the infinitesimal scale transformation Rt+ R+ dR, 
that is 


c= £0,c+ yOyc+ z0z€ (121) 
and isolating 0,c, yields 
Onc _ cos a + cos 3 (122) 
Cc x cos 8 cosy + ycosacosy + zcos acos B 


which, using al-Kashi’s theorem for the cosines and work- 
ing things out a little bit, can be reexpressed as 


Oz 
C G+ry (123) 
c (cty+z)\(a+ty—-z) 
and integrated out to give the general solution 
zr+y—zZ 
c(z,y, 2) = k(a, a 124 
(e.y.2)=Meyyfore = (124) 


where k(z,y) is a symmetric function of « and y. In 
the particular case of a right triangle with hypotenuse 
of length z, c = V2r where r is the inradius length. 
Computing the area of this right triangle yields r(x + 
y +z) = xy. Substituting r in the expression of c and 
inserting the latter in eq. (124), we find k(a#,y) = //xy. 
Hence 


uty-Z 


c(z,y, 2) =4/axy (125) 


e+ytz 


that is, the distance between the incenter and the foot of 
the y = (a,y) angle bisector. Dividing c(z, y, z) in this 
equation by the full length of the bisector d(x, y, z) in eq. 
(48) gives the ratio (x + y)/(c + y + z), meaning that 
the relative position/height of the incenter on the angle 
bisector (measured from its foot) is equal to the ratio 
of its corresponding sidelength and the perimeter of the 
triangle, that is z/(v+y+z); a result that can easily be 
verified in barycentric coordinates. 
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XVI. THE INRADIUS 


If we were able to do it for the circumradius, section [X] 
we should not doubt of Newton’s appetancy to do it for 
the inradius r! It would also be a smooth function of the 
sidelength x, y and z, that is 


(126) 


r=r(x,y,2) 


Strange as it may seem, since the position of the incen- 
ter matters, the circumcircle looks again to be the best 
place to do the job for the inradius. After a small dis- 
placement of the (x, y) vertex along the circumcircle, we 
have infinitesimal deviations dx and dy, while 6z = 0. At 
first order, 
or = Opr 6x + Oyr dy (127) 
dx and dy are determined in the same way as in the last 
sections — see for instance eq. for the case of the 
circumradius: 
dx = cosade 


— dy =cos Bde (128) 


For or we will refer to the geometrical considerations of 
section [XIV] and consider the small right triangle of legs 
of lengths —dr and 6s where ds is the distance travelled 
by the foot of the inradius perpendicular to the z-length 


side. First of all, we check that this inradius foot divides 
the z-length side in two segments of lengths (a — y + 
z)/2 and (—a + y+ z)/2 respectively, so that after the 
infinitesimal displacement it is moved of ds = (da—dy)/2. 
Secondly, in this small right triangle, the angle opposite 
to the —dr-length leg is equal to the angle between the 
internal y angle bisector and the height from the same 
vertex (or any line perpendicular to the opposite z-length 
side), that is y/2 — (7/2 — a) = (a — £)/2 since y = 
mt — (a+). We have 


— 6 da—6 
—6r = tan 5 * Ss 2 
= sina — sin ., (129) 
2 
thanks to eq. (128) and the development of the tangent 
a-6 sina-sinG 
a 2 cos a + cos 8 at30) 


Plugging in results (128) and into equation (127) 


gives a partial differential equation whose analogue can 
be obtained by considering a rotation of the (z, x) vertex 
instead of the (x, y) one. Both equations read 


sin 8 — sina 
2 
sina — siny 


2 


= cos a0,r — cos BOyr 131) 
= cos yOzT — COS AOzr 132) 
Combining those with the differential equation arising 


from the infinitesimal scale transformation Rw» R+ dR, 
that is 


r=«O0,r+yO0yr + 2027 133) 
and isolating 0,7, yields 


cos 3 cos y 
r 
x cos B cosy + ycosacosy + zcosacos G 
_ y(sin 6 — sina) cosy + 2(sin y — sin a) cos B 


= 134 
2(x cos B cosy + ycosacosy + z cos a cos 3) ss) 


Or? 


which, using the fundamental law for the sines, al-Kashi’s 
theorem for the cosines and working things out a little 
bit, can be reexpressed as 
F (x? — y? +. 2?) (a? + y? — 2?) 
r r 
x x(2(y2z? 4 272 + xy?) — 74 — y! = 24) 
_G-a)(@* +y? —27)4+(@— 2)(2? - y? +27) 
Qa /2(y2z? + x22? + wy?) — at — yt — 4 


(135) 


This equation is linear. Its general solution reads 


—0* + (y+ 2)a? = (y + 2)(y — 2)? + wey, 2) 


r(x, y,2) = 


2,/2(y222 f 2 72 f x7y?) 4 y* ZA 


where c(y, z) is a homogeneous function of y and z. Since 
r(x,y,z) must be symmetric in z, y and z, we have 
c(y,z) = y? + 27 + myz where m is a real constant that 
can easily be determined by considering, for instance, an 


13 


equilateral triangle with = y = z = 2/3r. We find 
m = —2. Hence, after factorization and simplification, 


_ f(-etyte)(e—yt+2z)\e+y—z) 
reewns)= A(a + y +2) 


(136) 
that is, the expected expression of the inradius as the 
area (see section [V) divided by twice the perimeter. 


XVII. THE ANGLE BISECTORS PROBLEM 


There is a well-known problem about the internal angle 
bisectors of a triangle: given three real and positive num- 
bers a, b, c, find the sidelength x, y, z of the triangle that 
would admit those numbers as angle bisector lengths. A 
lot of things have been said about it. We would not be 
far from the truth if we say that there is a solution, but 
that it is the solution of a polynomial equation of too 
high a degree to be solved [9]. See {I7}{19] for a histor- 
ical perspective and more recent attempt. We will just 
propulse Newton as a new challenger for a similar — but 
not equivalent — problem, that is identifying a, b, c not 
as the full lengths of the angle bisectors, but of their part 
going from the vertex to the incenter. This problem has 
a solution, but as it is the solution of a third degree poly- 
nomial equation, we will not try to write it down. Let 
Newton state it this way: each of the sidelengths must be 
a smooth function of a, b, c, the length of the segments 
which bisect the angles a, 3 and ¥y respectively, opposite 
to the sides of lengths x, y, z respectively. For z, we have 


z = 2(a, b,c) (137) 


Again, since the position and the movement of the incen- 
ter matter, the circumcircle is the best place to operate. 
After a small displacement of the (x,y) vertex along the 
circumcircle, we have the infinitesimal deviations da, 6b 
and dc, while 6z = 0. At first order, 

Oaz 6a + Opz 6b + O-z dc = 0 (138) 
If dp is the distance between the incenters before and 
after the small displacement, the angle between this dp- 
length segment and the —da-length is an inscribed angle 
of the circle along which the incenter travels, and has 


magnitude 6/2. Analogously, the angle between the dp- 
length segment and the db-length has magnitude a/2. 


- B 
da= cos 5 Op 


6b = cos 5 op (139) 
It is a bit more complicated for dc. We first observe that 
—dc is the length of one of the legs of a small right triangle 
with hypotenuse of length 62, the distance travelled by 
the moved vertex along the cicrumcircle. The magnitude 
of the angle between those two small sides is given by 
y/2+ 6. Since y = 7 — (a+ 8), the magnitude of the 
complemetary of this angle is (a — @)/2. Hence 


ap 
2 


dc = —sin be (140) 


To connect 6 to dp, just consider 6@ as the product of 
the circumradius length R by the corresponding central 
angle, that is twice the inscribed angle 6a = —d6 (since 
dy = 0). This inscribed angle da is itself connected to 
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dp by trigonometry: ada/2 is the second leglength of the 
small right triangle of first leglength —da and hypotenuse 
dp. But the circumradius length is not supposed to ap- 
pear in our equations; we can get rid of it in favour of z by 
reminding that siny = z/(2R), as can be seen by moving 
the (x,y) vertex along the circumcircle, thus preserving 
y, until the x or y-length side coincides with a diameter. 
These considerations can be summarized as follows: 


) 
=tRiy ¢° =eHe oy oR= Gal) 
2 2 sin 7 
Inserting those equations in eq. (140) yields 
- 2 
dc = —sin* ; e — sin ; Sp (142) 


We need a little trigonometry to move forward. Apply- 
ing the law of sines in at least two of the three triangles 
generated by the a, b and c-length segments in the initial 
triangle, we find 

B 


;, @ , . 
asin ~ = bsin= = esin — 
2 2 


2 


Furthermore, if we apply the al-Kashi theorem to the 
angle (a,b) =a — (a+ B)/2 =7/2+4+ 7/2, we have 


(143) 


y ga? —P 


i = 144 
sas) 2ab ve 
Inserting this in eq. (143) gives 
a c22—-a—b B c2-ad-—b 
5 = i = 145 
2 aad Moy pay 


——_. 


Applying the al-Kashi theorem to the angles (a, z) = a/2 
and (b, z) = 8/2 respectively, we find 


a 2+a?—6? B 2-d@s+b 
Cos = cos — 


2 2az 2 2bz 
Inserting those results in the development of cos(y/2) = 


sin(a/2 + 6/2) = sin(a/2) cos(8/2) + cos(a/2) sin(6/2) 
yields 


(146) 


y cz 227-a?—0? 
2 ab 2ab 


Note that submitting sin(7/2) in eq. and cos(y/2) 
in eq. to the fundamental law leads to a polynomial 
equation of degree 3 in z?, that is, the complicated but 
existing solution to our problem. 

But in this quest, we were following Newton on another 
path: it is now time to exploit the expressions of the small 
deviations in eq. and — developing sin(a/2— 
8/2) and using sin y = 2 sin(y/2) cos(y/2) in the latter — 
and plug in them into eq. to finally obtain our first 
partial differential equation 

za? +b? 4 2+a?— 0? b? — a? 


Th Be 2a ere 


(147) 


0.z =0 

(148) 
Combining it with the differential equation arising from 
the infinitesimal scale transformation R'> R+ dR, that 
is 


Z=A0,2 + bOnz + 60.2 (149) 


and isolating 0,2 and Oz, yields 


e2tat—b? ¢ 22-a?4+b? 


On% = OeZ (150) 
2az Zz 2az 
eat 6 24+a2— 
= P 151 
Ovz 2bz z 2bz Oc2 et) 


That was the easy part. We need a third partial dif- 
ferential equation to complete the system and be able 
to solve it. The analogue of eq. (148) can be obtained 
by considering a small displacement of the (y, z) vertex 
along the circumcircle, instead of the (x,y) one, but it is 
more complicated since now 6z # 0. At first order 


62 = Ogz 6a + Oz 6b + 0,2 0€ (152) 
Following a similar reasoning as previously we find 
lt a eae ne 
= 1 
da sin — mat sin 5 Op 53) 
6b = —cos ; op 154) 
oe ee 
dc = cos . op 155) 
2 
bz = —cosy aa 5 sin 3 op 156) 
where the latter derives from dz = — cosy 6, by analogy 


with what was done, for instance, to determine dy in eq. 
(66). Using the same trigonometry as for the previous 
case, plus cosy = cos?(7/2) — sin?(7/2), and plugging in 
the results into eq. (152), we have 


2,2 272 »2_ g2 _ p2 


cz" —a@ 
ab2c 2a 
cz 22 —- a — ab 2 -— a? +0? 
= Ogz 
ab 2b CZ 2b 
ee ee © ane ane 
ane Onz —— On% (157) 
ab 2a 2z 


All that remains is to replace, in this equation, 0,z and 


Oz by their expression in eq. (150) and (151) respec- 
tively. After a tedious calculation we finally obtain the 
differential equation 


COZ = (158) 
2(z? — a? — b*)(—24 + 2(a? + b2)z? — (a? — b?)?) 


76 — 3(a? + b?) 24 + 3(a? _ b?)? 2? = (a? +4 b) (a? _ b?)? 


which can be integrated out to give the general solution 


z?)(z? — (a— b)?) 
2(z2 — a? — b?) 


=k(a,b)c (159) 


where k(a,b) is a symmetric function of a and b. Here 
z(a,b,c) si implicitely given as a root of a polynomial 
equation of degree 3 in z?... And this equation happens 
to be the same as the one we have mentionned earlier if 
we take k(a,b) = £1/(ab). 
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XVIII. THE POLYNOMIALS 


One of the lessons that we can draw from this last at- 
tempt is that the method cannot evade the sometimes 
complicated question of finding the roots of a polynomial. 
Newton surely was interested in polynomials. Could he 
have tried to find their roots through calculus? Let us 
first look at a polynomial P(x) of degree 2 with real co- 
efficients a, b and c. Each of its roots, if any, satisfies the 
equation 

ax? +ba+ce=0 (160) 


Assume that the root x of our polynomial is a smooth 
function of the coefficients: 
x = x(a, b,c) (161) 


After an infinitesimal deviation of those coefficients, we 
have, at first order: 
6x = Ogx da + Opx 6b + O,x bc (162) 


But eq. 
order, 


(160) must also remain valid, so that, at first 


xa +26b + dc = —(2ax + b)da (163) 


If we consider first dc £ 0 with da = 6b = 0, then 6b 4 0 
with da = 6c = 0, and finally 6a 4 0 with 6b = dc = 0, 
we have repsectively 


2ax +b 
2 


2ax +b 


dc = —(2ax 4+ b)da 6b = — dx da=— bx 
that we plug in into eq. (162) to obtain three partial 


differential equations 


—1 —2x —x? 


2ax +b one 2ax +b = 2ax +b en) 


cl 


The general solutions of these equations are respectively 


ax* + br +c=k,(a, b) (165) 
(2ax + b)? = b? + ky(a,c) (166) 
(2ax + b)* — b? = ka(b, c)a (167) 


where k.(a,b), ky(a,c), and ka(b,c) are arbitrary func- 
tions. Using the two last solutions leads to ky(a,c) = 
ka(b,c)a = k(c)a where k(c) is an arbitrary function; 
comparing this to the first solution yields 4k,(a,b) = K 
and k(c) = K — 4c, with K a constant that must vanish 
since x is a root of the polynomial; which brings us back 
to eq. (160). 

For a polynomial of degree n with real coefficients, the 
roots satisfy the polynomial equation 


(168) 


Assume that the root z is a smooth function of the coef 
ficients: 


(169) 


X= U(do,.--, An) 


After an infinitesimal deviation of those coefficients, we 
have, at first order: 


n 


6x = S- Oa, t Ok 


k=0 


(170) 


If we replace a, by a, +daz and x by x+6z in eq. (168), 
we have at first order again 


n 


n 
y x bay = — y iaja' | dx 
k=0 


i=0 


(171) 


Now, by successively keeping only one non-zero coeffi- 
cient deviation da; for k = 0,...,n, we find the (n + 1) 
relations 


n 
x ba, = — ) ia; | bx 


i=0 


(172) 


that can be plugged in into eq. (170) to yield the (n+ 1) 
partial differential equations 


(>: ie) Oa, & + xc =0 


(173) 
1=0 


for k = 0,...,n. The general solution of this system 
should provide the roots of the polynomial! Each of these 
equations can be expressed as an ordinary first-order dif- 
ferential equation that turns out to be exact. For each 
k, we can indeed see our polynomial as a function ux of 
x and ag, i.e. 


Ug (2, ap) = S- ajax" (174) 
i=0 
Its partial derivatives 
OvUR = \ ae and 0g, Un = x* (175) 


i=0 


are naturally the coefficients of eq. (173) so that for each 
k, we have the total differential du, = 0, and ug(a, ag) 
must be equal to an arbitrary function of all the a; for 
i #k. Equating those (n + 1) arbitrary functions forces 
them to be constant and eventually to vanish since x is 
a root of the polynomial. We end up with eq. (168). 
In other words, we have found a set of partial differen- 
tial equations equivalent to the polynomial equation, but 
the solution of this set is naturally given in the implicit 
form of the polynomial equation itself, as we could have 
expected! Newton, at least, would have... 


XIX. ARCHIMEDES OF SYRACUSE 


We don’t need to speculate on Newton’s aspirations, to 
pretend to discover that the method can be used to com- 
pute the area of a circle and the volume of a sphere as 
functions of their radius length r by observing the way 
they behave under the transformation r+ r + dr. 
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For the circle area A(r), we have 6A = A’dr, and it 
is easy to calculate 6A = 2zrér, so that A’ = 27r and 
A(r) = mr? +k where k is a constant equal to 0 since 
A(0) = 0. 

We cannot do it for the perimeter P(r) of the circle 
since the intermediate result 6P = 276r that leads to 
P’ = 27 and P(r) = 2z7, precisely relies on the fact that 
P = 2nr, that is, the definition of 7. 

If V(r) is the volume of the sphere, then dV, = Vjor. 
By imparting a small thickness dr to the spherical en- 
velope of area 4ar?, we calculate 6V; = 4rr?6r, so that 
Vi = 4rr? and V,(r) = 4rr?/3+k where k is a constant 
equal to 0 since V,(0) = 0. 

We could do it for the area of the sphere A,(r), but 
we did not find a way to compute dA, = 8mrér with- 
out using integral calculus and/or, at least, the cylindri- 
cal projection that Archimedes used to directly compute 


A,(r) = 4nr? PO]. 


XX. A SIMPLE EQUATION 


The equations (70), (81), (94), (121), and 


that were used to determine the circumradius, the law of 
sines, etc. could be formulated in a more general way — 
still specific to Euclidean Geometry: 


iv 
J Liste) = me MeiOof (176) 


where f is a metric quantity (length, area, volume,...) 
and n its dimension (1,2,3,...), while x; are the metric 
quantities on which f depends, and n,; their own dimen- 
sions. Every theorem of Euclidean Geometry expressed 
as a smooth function, every length, surface, volume, hy- 
persurface formula should satisfy equation which 
insures that all metric quantities involved in it are dimen- 
sionally correctly mixed together. It is naturally true for 
all the theorems re-derived in this article. 

Besides the six cases mentionned above, we show how 
we could have used equation to derive Pythagoras 
and Heron’s theorems. In the first subsection we observe 
how this equation specifically acts by restricting a large 
class of solutions. The second subsection is devoted to 
the study of a n = 2 case, since Heron deals with the 
area of triangles. 


A. Pythagoras 


Imagine that instead of linearily increasing x (and then y) 
of an infinitesimal 62, we would have chosen to perform 
an infinitesimal rotation of the hypotenuse around, say, 
the (y, z) vertex. Here 6z = 0, we would then have 

0,2 6% + Oyz dy = 0 (177) 
where the small deviations dz and —dy happen to be the 
leglengths of a right triangle similar to the initial one — 
since its hypotenuse is, at first order, perpendicular to 


the hypotenuse of the initial triangle —, yielding 


by = da (178) 
Yy 
y y 
—dy 
x ox x 
Plugging in this equation into equation (177) gives 
YOu — LOyz = 0 (179) 


which admits a large class of solutions of the form 
z(x,y) = f(a? + y”). To complete the derivation of the 
Pythagorean theorem, we need to consider the way z be- 
haves, at first order, under a scale transformation, that is, 
Rt R+O6R, where R is the circumradius (here R = z/2) 
or any other length to be used as a scale. We have 


Z=LOpz + Y Oyz (180) 


which corresponds to equation (176). Joining equations 
(179) and (180) in a system and isolating 0,z yields 


Onz x 
eS ys 181 
z x? + y? eh) 
and its general solution 
2(r,y) = k(y) Va? + y? (182) 


with k(y) a function of y. But the equation (181) has its 
analogous for 0,z, generating by integration a function 
k(x) which must be equal to k(y), that is, to a real and 
positive constant k. Considering a triangle degenerated 
in a segment implies that k = 1. 


B. Heron 


We still postulate that 


A= A(a,y,2) (183) 
but we now consider the triangle in its circumcircle. If 
the (x, y) vertex is slightly moved along the circumcircle, 
it generates the infinitesimal deviations dx and dy, while 
dz = 0. At first order, we have: 

6A = 0,A dx + OyA dy (184) 
The deviations 62 and dy are the same as in the cir- 
cumcircle section, equation while 6A = zdh/2. To 
determine 6h, we observe that it is one of the leglengths 
of a small right triangle of hypotenuse of length 62, and 
that this right triangle is similar to another one of hy- 
potenuse of length R and leg with corresponding length 
d,, the distance between the circumcenter and the height 
of length h, that is, the distance between the middle of 
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the z-length side and the foot of the h-length height, or 
else 


; (185) 


This can be seen by sub- 
tracting t, in eq. (29), to 
z/2. Hence 

d, 7 y? — x? 
— 2R 


5 (186) 


6h = = 
R 

_yp-2 

~ AR 


oe 


>06A 


But in the circumcircle, 
the angle between the h- 
length height and the z- 
length side is the same as 
the one between the y- 
length side and the circum- 
radius from de (x,y) ver- 
tex, thus h/a = y/(2R). 


(187) 
Thus 


be (188) 
Plugging in equations and (188) into equation (184), 
simplifying by 6@/z and re-arranging the terms leads to 


y2 —a2 +2? 
y2 — x2 


a2 —y? +2? 


2A=2 ae 


OA 


d,A (189) 


This equation can be joined to its symmetric couterpart 
corresponding to the infinitesimal displacement of the 
(z,x) vertex along the circumsircle, while dy = 0, that is 
a2 — 224 y? 22 

pei z o2 — 


2A=z 


+y? 
2 O,A (190) 
and the equation originating from the infinitesimal scale 
transformation RH R+dR 
2A=20,A+yOA+z0,A (191) 
in a system of three equations that we can solve in 0,A 
to obtain the partial differential equation 
O,A 1 —4a? + 4a(y? + 27) 
A 2 [—a4 — yt — 24 + 2(y?2? 


x222 + x2y2)| 
(192) 
which can be integrated out to give the general solution 
A(x, y, 2) (193) 

= by, 2) =a FF GPE aE a) 
where k(y, z) is an function of y and z. Since the function 
A must be symmetric in x, y and z, k(y, z) = c where cis 
a constant. We just have to compute the area of, say, an 
isocele right triangle of leglength 1, that is A(1,1, /2) = 


1/2, to find out that c = 1/4. After factorization, the 
result is equation (34), that is, Heron’s theorem. 


CONCLUSION 


With this paper, we present an alternative way to de- 
rive classical theorems in Euclidean geometry. Not all 
theorems, of course. It does not work for theorems in 
discrete geometry or involving number theory, for theo- 
rems stating that this or that line cuts another at this 
or that point, is perpendicular or tangent to this or that 
circle, etc. It has to be a theorem involving an equation 
that defines a function, which will be seen as a particu- 
lar solution of a (system of) differential equation(s). The 
proofs that we propose are not necessarily simpler than 
others. They do not evade the geometric difficulties at 
stake. We displace the argument of the proof into the 
game of infinitesimals, but it remains as geometric. 

The main advantage of this method is that the theorem 
does not need to be known. We start with a function, 
any function, and observe the way it behaves, at first 
order, under small deviations of some quantities. In the 
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best case, it gives us a (system of) differential equation(s) 
that we can solve and, therefore, discover the theorem. 
We chose to use it to rediscover about 20 theorems or 
identitites in the long history of Euclidean geometry, but 
we hope it may be used to discover new theorems, per- 
haps in other fields of mathematics. 
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